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Abstract 

We compute the next-to-leading order QCD corrections to spectator effects in the 
hfetime ratios of beauty hadrons. With respect to previous calculations, we take into 
account the non vanishing value of the charm quark mass. We obtain the predictions 
T{B+)/T{Bd) = 1.06±0.02, T{Bs)/T{Bd) = l.OOiO.Ol and r(Afe)/T(Sd) = 0.90±0.05, 
in good agreement with the experimental results. In the case of t{Bs) /T[Bd) and 
t[Ki))/t{B(i), however, some contributions, which either vanish in the vacuum inser- 
tion approximation or represent a pure NLO corrections, have not been determined 
yet. 



1 Introduction 



The inclusive decay rates of beauty hadrons can be computed by expanding the amphtudes 
in increasing powers of kqcDl'^h [0, H] and using the assumption of quark-hadron duahty. 
The leading term in this expansion reproduces the predictions of the naive quark spectator 
model, and the first correction of (9(l/mb) is absent in this case. 

Within this theoretical framework, up to terms of 0{l/ml) only the b quark enters 
the short-distance weak decay, while the light quarks in the hadron interact through soft 
gluons only. For this reason, the lifetime ratios of beauty hadrons are predicted to be unity 
at this order, with corrections which are at most of few percent. 

Spectator contributions, which are expected to be mainly responsible for the lifetime 
differences of beauty hadrons, only appear at 0{l/ml) in the Heavy Quark Expansion 
(HQE). These effects, although suppressed by powers of l/rrib, are enhanced, with respect 
to leading contributions, by a phase-space factor of IGvr^, being 2^2 processes instead of 
1 — > 3 decays. Indeed, the inclusion of these corrections at the leading order (LO) in QCD 
has allowed to reproduce the observed pattern of the lifetime ratios p, ^ . 

In this paper we compute the next-to-leading order (NLO) QCD corrections to spectator 
effects in the lifetime ratios of beauty hadrons. In the limit of vanishing charm quark mass, 
these corrections have been already computed in ref. 0, and we refer to this paper for many 
details of the NLO calculation. We find that the inclusion of the NLO corrections improves 
the agreement with the experimental measurements. Moreover, it increases the accuracy 
of the theoretical predictions by reducing significantly the dependence on the operator 
renormalization scale. 

After our calculation was completed, the NLO corrections to spectator effects in the 
ratio T{B^)/T{Bd) have been also presented in ref. We agree with their results. With 
respect to ref. 0, we also perform in this paper the NLO calculation of the Wilson coeffi- 
cients entering the HQE of the ratios t{Bs)/t{Bci) and r(Afe)/r(i?rf). 

By using the lattice determinations of the relevant hadronic matrix elements [[f|-||10 



we have performed a theoretical estimate of the lifetime ratios and obtained the NLO 
predictions 

r{B+) t(B,) T(Ab) 

— f = 1.06 ±0.02, = 1.00 ±0.01, = 0.90 ± 0.05 . 1 

r{Bd) T{Bd) T{Bd) ^ ' 



These estimates must be compared with the experimental measurements [p 



= 1-074 ±0.014, 4^ = 0.948 ± 0.038 , = 0.796 ± 0.052 . (2) 

T{Bd) T{Bd) T{Bd) 

The theoretical predictions turn out to be in good agreement with the experimental data. 
As we will discuss below, however, in the case of t{Bs) /T{Bd) and T{Ab)/T{Bd) the theo- 
retical predictions can be further improved, since some contributions, which either vanish 
in the vacuum insertion approximation (VIA) or represent a pure NLO correction, are still 
lacking. With respect to the results of ref. we find that the NLO charm quark mass 
corrections are rather large for some of the Wilson coefficients. Nevertheless, the total 
effect of these contributions on the lifetime ratios is numerically small. 
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An important check of the perturbative calculation performed in this paper is provided 
by the cancellation of the infrared (IR) divergences in the expressions of the Wilson co- 
efficients, in spite these divergences appear in the individual amplitudes. The presence 
of these divergences explicitly shows that the Bloch-Nordsieck theorem does not apply in 
non-abelian gauge theories llll-ll^. We have also checked that our results are explicitly 
gauge invariant and have the correct ultraviolet (UV) renormalization-scale dependence as 
predicted by the known LO anomalous dimensions of the relevant operators. 

The HQE of the lifetime ratios results in a series of local operators of increasing dimen- 
sion defined in the Heavy Quark Effective Theory (HQET). Indeed, renormalized operators 
in QCD mix with operators of lower dimension, with coefficients proportional to powers 
of the 6-quark mass. In this case, the dimensional ordering of the HQE would be lost. In 
order to implement the expansion, the matrix elements of the local operators should be 
cut-off at a scale smaller than the 6-quark mass, which is naturally realized in the HQET.Q 
In the case of t{B^)/t{B(i), only flavour non-singlet operators enter the expansion. For 
these operators, the mixing with lower dimensional operators is absent, and the HQE can 
be expressed in terms of operators deflned in QCD. For these operators, the matching 
between QCD and HQET has been computed, at the NLO, in ref. 

As mentioned before, the NLO predictions for the lifetime ratios T{Bs)/T{Ba) and 
r(Ab)/r(i?rf) may still be improved in some ways. In particular: 

- the contributions of the four-fermion operators containing the charm quark fleld and 
of the penguin operator (see eqs. (|r^) and ([TTp) have not been determined yet. These 
contributions, which either vanish in the VIA (charm operators) or represent a pure 
NLO correction (penguin operator), do not affect the theoretical determination of 
t[B^) /r^Bd) and represent only an S'f/(3)-breaking effect for t{Bs) /t^Bcl)] 

- the lattice determination of the hadronic matrix elements have been performed by 
neglecting penguin contractions (i.e. eye diagrams), which exist in the case of flavour 
singlet operators. Also in this case, the corresponding contributions cancel in the 
ratio T{B^)/r{Bd) but affect r{Ah)/T{Bd) and t{Bs) /t{B(i), the latter only through 
SU (3)-breaking effects; 

- the NLO anomalous dimension of the four-fermion DB = operators in the HQET 
is still unknown. For this reason, the renormalization scale evolution of the A;, matrix 
elements, computed on the lattice in the HQET at a scale smaller than the 6-quark 
mass p, y], has been only performed at the LO. For B mesons, a complete NLO 
evolution has been performed by using operators deflned in QCD ||10|| . 

For these reasons, at present, the best theoretical accuracy is achieved in the determination 
of the lifetime ratio T{B^)/T{Bd). 

We conclude this section by presenting the plan of this paper. In sect, ^we review the 
basic formalism of the HQE applied to the lifetime ratios of beauty hadrons. Details of the 
NLO calculation and the numerical results obtained for the Wilson coefficients are given 
in sect. ^. In sect. ^ we present the predictions for the lifetime ratios of beauty hadrons. 

^On this point, we disagree with the general statement of ref. |^ according to which the HQE can be 
performed in QCD also in the presence of mixing with lower dimensional operators. 
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Finally, we collect in the appendix the analytical expressions for the Wilson coefficient 
functions. 



2 HQE for the lifetime ratios of beauty hadrons 

Using the optical theorem, the inclusive decay width T{Hf,) of a hadron containing a b 
quark can be written as 

m,) = j^lm{H,\r\H,) , (3) 



where the transition operator T is given by 

r = i I d'^xT{n 



eff ' 



(0) 



(4) 



and 'H^j^^^ is the effective weak hamiltonian which describes AB = 1 transitions. 

By neglecting the Cabibbo suppressed contribution oib ^ u transitions (|V^6p/|V^6p ~ 
0.16A^) and terms proportional to |V^(i|/|V^s| in the penguin sector, the AB = 1 effective 
hamiltonian can be written in the form 



o/AB=l 
^eff 



j=3 l=e,fj,,T 

The Ci are the Wilson coefficients, known at the NLO in perturbation theory 
the operators Qi are defined as 



+ h.c.. (5) 

[0, and 



Ql = {biCj)v~A{Ujdi)v-A, 
Ql = {biCj)v-AiCjSl)v-A, 

Qs = {biSi)v-A'^{qjqj)v-A , 

Q5 = (biSi)v-Aj2(Wj)v+A , 



Q2 = {biCi)v~A{Ujdj)v-A, 
Q2 = ibiCi)v-A{CjS-)v-A, 

Qa = {biSj)v-A'^{qjqi)v-A, 
q 

Qg = ibiSj)v-A^{qjqi)v+A , 



(6) 



Ql = ibiCi)v-AM 



V-A 



where d' = cos 9 cd + smdcS and s' = — sin 6cd + cos 6 cS {6c is the Cabibbo angle). Here 
and in the following we use the notation {qq)v±A = ^7^(1 ^l5)q and (qq)s±p = q{l ±75)^. 
A sum over repeated colour indices is always understood. 

Because of the large mass of the b quark, it is possible to construct an Operator Product 
Expansion (OPE) for the transition operator T of eq. (^, which results in a sum of local 
operators of increasing dimension |^. We include in this expansion terms up to 0{l/ml) 
plus those l/mf corrections that come from spectator effects and are enhanced by the phase 
space. The resulting expression for the inclusive width of eq. @ is given by 



1927r3 



.(3) 



2M 



+ c 



(5) 



gs {ba^uG^''b)H, , 967r 



2M. 



+ 



mi 



E 



.(6) 



2M, 



Hb 



(7) 
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Figure 1: Examples of LO contributions to the transition operator T (left) and to the 
corresponding local operator (right). The crossed circles represent the insertions of the 
AB = 1 effective hamiltonian. The black squares represent the insertion of a AB = 
operator. 

where (■ ■ ■)h^ denotes the forward matrix element between two hadronic states Hh, defined 
with the covariant normahzation. The operators O^^"* are a set of four-fermion, dimension- 
six operators to be specified below, which represent the contribution of hard spectator 
effects. At the lowest order in QCD, the diagrams entering the calculation of T{Hh) are 
shown in fig. ^ 

The matrix elements of dimension-three and dimension-five operators, appearing in 
eq. (1^, can be expanded by using the HQET 

= 2MH.(l- ^-™;,^^'^^> +0(l/n.a). 

J.(Sff„.G""'())„. = 2M„, (2i^UK) + 0(l/m(,)) . (8) 

By using these expansions, we can finally compute, from eq. (0), the lifetime ratio of two 
beauty hadrons 

T{Hl,) 2ml V ) ml 
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Z^^'^ on#.. on#... • 



b 



The dimension- six operators in eq. which express the hard spectator contributions, 
are the four current- current operator^ 



(10) 



OI = {bq)v-A iqb)v-A , 0| = {bq)s-P ms+P , 

OI = ibt^q)v-A {qt%)v-A , OI = {bt-q)s-p iqt%)s+p , 
with q = u,d, s, c, and the penguin operator 

Op = {bt%)v E i.(lt»v- (11) 



q=u,d,s,c 



In these definitions, the symbols b and b denote the heavy quark fields in the HQET. 

The Wilson coefficients c^^^ and c^^-* , of the dimension-three and dimension- five operators 
in eq. (H), have been computed at the LO in ref. [jl8|, while the NLO corrections to c*-^-* have 



been evaluated in The NLO corrections to c^^-* are still missing. Their numerical 

contribution to the lifetime ratios, however, is expected to be negligible. 

The coefficient functions of the dimension-six current-current operators Ol, have been 



computed at the LO in ref. ^ for q = u,d, s, and in ref. p6| for q = c. The coefficient 
function of the penguin operator Op vanishes at the LO. 

In this paper we have computed the NLO QCD corrections to the coefficient func- 
tions of the operators Ol with q = u,d, s. The operators containing the charm quark 
fields contribute, as valence operators, only to the inclusive decay rate of Be mesons, and 
their contribution to non-charmed hadron decay rates is expected to be negligible. The 
calculation of the NLO corrections to these coefficient functions, as well as the NLO cal- 
culation of the coefficient function of the penguin operator, has not been performed yet. 
The non-perturbative determinations of the corresponding matrix elements are also lack- 
ing at present. However, as we will discuss in sect. ^, these contributions only enter the 
theoretical estimates of the ratios t{Bs) /t{Bii) and r(Afe)/r(i?^), and vanish in the VIA. 



3 NLO calculation of the Wilson coefficients 

All the details of the matching procedure used to determine the Wilson coefficients at 
the NLO have been given in ref. [^, where the calculation has been performed in the 
limit of vanishing charm quark mass. A finite value of the charm quark mass does not 
introduce conceptual difficulties in the calculation, besides requiring the evaluation of one- 
and two-loop integrals with an additional mass scale. For this reason, we only remind in 
this section the general strategy of the perturbative calculation, and refer the interested 
reader to ref. for further details. In this section, we will also present the numerical 
results for the Wilson coefficients, while the analytical expressions of these coefficients are 
collected in the appendix. 
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Figure 2: Feynman diagrams which contribute at NLO to the matrix element of the tran- 
sition operator T in the case q = s. In the other cases, q = u, d, diagrams D14 and D15 
are Cabibbo suppressed and have been neglected in the calculation. 
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Figure 3: Feynman diagrams which contribute, at NLO, to the matrix element of the AB = 
operators entering the HQET. 



In order to compute the Wilson coefficients of the AB = operators at the NLO, we 
have evaluated in QCD the imaginary part of the diagrams shown in fig. Q {full theory) and 
in the HQET the diagrams shown in fig. ^ {effective theory). The external quark states 
have been taken on-shell and all quark masses, except ruh and rric, have been neglected. 
More specifically, we have chosen the heavy quark momenta pi = ml in QCD and fc;, = in 
the HQET, and Pq = for the external light quarks. In this way, we automatically retain 
the leading term in the l/rub expansion. We have performed the calculation in a generic 
covariant gauge, in order to check the gauge independence of the final results. Two-loop 
integrals have been reduced to a set of independent master integrals by using the recurrence 
relation technique [^-[^ implemented in the TARCER package [^. Equations of motion 
have been used to reduce the number of independent operators. 

Some diagrams, both in the full and in the effective theory, are plagued by IR diver- 
gences. These divergences do not cancel in the final partonic amplitudes, and provide an 
example of violation of the Bloch-Nordsieck theorem in non-abelian gauge theories fT^ - 
IR poles, however, are expected to cancel in the matching and we have explicitly 
checked that, in the computation of the coefficient functions, this cancellation takes place. 

We use Z)-dimensional regularization with anticommuting 75 (NDR) to regularize both 
UV and IR divergences. As discussed in details in ref. 0, the presence of dimensionally- 
regularized IR divergences introduces subtleties in the matching procedure. The matching 
must be consistently performed in D dimensions. This requires, in particular, enlarging 



^Note that, with respect to ref. we use a different basis of operators. 
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the operator basis to include (renormalized) evanescent operators, which must be inserted 
in the one-loop diagrams of the effective theory. Because of the IR divergences, the matrix 
elements of the renormalized evanescent operators do not vanish in the D — > 4 limit pl| , 
and give a finite contribution in the matching procedure. 

As a check of the perturbative calculation, we have verified that our results for the 
Wilson coefficients satisfy the following requirements: 

• gauge invariance: the coefficient functions in the MS scheme are explicitly gauge- 
invariant. The same is true for the full and the effective amplitudes separately; 

• renormalization- scale dependence: the coefficient functions have the correct logarith- 
mic scale dependence as predicted by the LO anomalous dimensions of the AB = 
and AB = 1 operators; 

• IR divergences: the coefficient functions are infrared finite. We verified that the 
cancellation of IR divergences also takes place, separately in the full and the effective 
amplitudes, for the abelian combination of diagrams . 

The AB = effective theory is derived from the double insertion of the AB = 1 effective 
hamiltonian. Therefore, the coefficient functions cl of the AB = effective theory depend 
quadratically on the coefficient functions Ci of the AB = 1 effective hamiltonian, and we 
can write 

4(/^o) = EC.(/^i)C,(/xi)F,^,//xi,/xo). (12) 



Eq. (jT^) shows explicitly the scale dependence of the several terms. We denote with /ii the 
renormalization scale of the AB = 1 effective hamiltonian, whereas Hq is the renormaliza- 
tion scale of the AB = operators. The coefficients Fj}- depend on the renormalization 
scheme and scale of both the AB = and AB = 1 operators. The dependence on the scale 
fii and on the renormalization scheme of the AB = 1 operators actually cancels, order by 
order in perturbation theory, against the corresponding dependence of the AB = 1 Wilson 
coefficients Ci. Therefore, the coefficient functions cl only depend on the renormalization 
scheme and scale of the AB = operators. We have chosen to renormalize these operators. 



in the HQET, in the NDR-MS scheme defined in details in ref. ^2 . 

The four fermion operators, in the effective AB = 1 hamiltonian, are naturally ex- 
pressed in terms of the weak eigenstates d' and s'. For this reason, we can write the 
coefficients Fj^^j of eq. ([I^ ) in the form 



F^^,^ = cos' + sin^ 9,F^;,^ 

F^,^ = sin' + cos^ e,F^;^^ (13) 

for i,j = 1,2. In the case q = s, the coefficient functions also receive contributions from 
the insertion of the penguin and chromomagnetic operators (diagram Di^ of fig. |]). Since 
the Wilson coefficients C^-Cq are small, contributions with a double insertion of penguin 



operators can be safely neglected. As suggested in a consistent way for implementing 
this approximation is to consider the coefficients C^-Cq as formally of 0{as)- Within this 
approximation, only single insertions of penguin operators need to be considered at the 
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q = d 


q = u 


q = s 




LO 


NLO 


NLO 


LO 


NLO 


NLO 


LO 


NLO 


NLO 






(me = 0) 






(m, = 0) 






(m, = 0) 






-0.02 


-0.03 


-0.03 


-0.06 


-0.33 


-0.29 


-0.01 


-0.03 


-0.03 


4 


0.02 


0.03 


0.03 


0.00 


-0.01 


-0.01 


0.02 


0.03 


0.04 


4 


-0.70 


-0.65 


-0.67 


2.11 


2.27 


2.34 


-0.61 


-0.51 


-0.56 


4 


0.79 


0.68 


0.68 


0.00 


-0.06 


-0.05 


0.77 


0.61 


0.64 



Table 1: Wilson coefficients c| computed at the LO and NLO, in the latter case with and 
without the inclusion of charm quark mass corrections at 0{as)- As reference values of the 
input parameters, we use = fii = = 4.8 GeV and rric = 1.4 GeV. 

NLO, and we can write 

4= E QC,F^^,^+2-^C2CsGPk,28 + 2 E J2C^CrPk,ir, (14) 
i,j=l,2 i=l,2r=3,6 

which generahzes eq. (|T^) for the case q = s. 

The analytical expressions of the coefficients F^ ^ - and Pfc,ij will be given in the appendix. 
For illustration we present here, in table |I], the numerical values of the coefficients c\ both 
at LO and NLO and, in the latter case, with and without the inclusion of the 0{as) 
charm quark mass corrections computed in this paper. As reference values of the input 
parameters, we use fio = fii = rub = 4.8 GeV and nic = 1.4 GeV. 

By looking at the results shown in table |l|, we see that the NLO charm quark mass 
corrections are rather large for some of the Wilson coefficients, although the total effect 
of these contributions on the lifetime ratios will be found to be small. The numerical ex- 
pressions of the lifetime ratios as a function of the 5-parameters will be given in the next 
section (eq.(p6D). In these expressions we will also include an estimate of the theoretical 
error on the coefficients, coming from the residual NNLO dependence on the renormal- 
ization scale Hi and from the uncertainties on the values of the charm and bottom quark 
masses and the other input parameters. 

The combinations c^ — cf (fc = 1, . . . 4) of Wilson coefficients, which enter the theoretical 
expression of the ratio T{B~^)/T{Bd), have been also computed, at the NLO, in ref. 0. In 
this case, the corresponding operators are the flavour non-singlet combinations (9^ — Of, 
which do not mix, with coefficients proportional to powers of the b quark mass, with 
operators of lower dimension. For this reason, the HQE can be also expressed in this case 
in terms of operators defined in QCD, and this is the choice followed by ref. 0. 
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In order to compare our results with those of ref. at the NLO, it is necessary to 
implement the matching, at 0{as), between QCD and HQET operators. This matching 
can be written in terms of the Wilson coefficients, in the form 



QCD 



1 + 



An 



S 



{nib) - cf{mb) 



HQET 



(15) 



where a common renormalization scale fi = rrib has been chosen for all the coefficients. 
The matrix S has been computed in ref. ||^ . It depends on the renormalization schemes of 
both QCD and HQET operators. In this paper, we have chosen to renormalize the HQET 
operators in the NDR-MS scheme of ref. |32]. By choosing for the QCD operators the 
NDR-MS scheme of ref. p , defined in ref. ||3^ , one finds that the matrix S is given by f\ 



S 



V 



32/3 





22/9 


1/3 \ 


-16/3 


-16/3 


4/9 


-2/9 


11 


3/2 


-7/2 


1/4 


2 


-1 


3 


5/2 ) 



(16) 



By using eq. (plSj), we have verified that our results for the combinations — cf agree 
with those of ref. p. Note that in the notation of |^ the labels u and d are interchanged 
with respect to our convention and that the Wilson coefficients are defined with a relative 
factor of 3. The numerical comparison of our results with those shown in table 1 of ref. IQ, 
however, shows some differences, particularly at the LO. The reason is that some of the 
coefficient functions, because of large cancellations, are extremely sensitive to the value of 
the coupling constant as{fi). 



4 Theoretical estimates of the hfetime ratios 

In this section we present the theoretical estimates of the lifetime ratios of beauty hadrons, 
obtained by using the NLO expressions of the Wilson coefficients and the lattice determi- 
nations of the relevant hadronic matrix elements []7|-||1CI||. 

The HQE for the ratio of inclusive widths of beauty hadrons is expressed by eq. (^, up 
to and including 1 / mf spectator effects. The combinations of hadronic parameters entering 
this formula at order 1/ml can be evaluated from the heavy hadron spectroscopy and 
one obtains the estimate 

^ = l-00-Ajl, ^ = 1.00-AJ;,, ^ = 0.98(1) -A^,,,, (17) 
where the As represent the 1 / ml contributions of hard spectator effects 

spec ^3c(3)Z. k I 2M^^ 2Mb, 



■^This matrix differs from the matrix C'l given in ref. ||^ by a Hnear transformation, since a basis of 
operators different from eq. (|o|) has been chosen in that paper. More specifically, S = —MC'[M~^, where 



M is defined in ref. 1^. 
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Note that, by neglecting these contributions, the predictions of eq. (|T^) are incom- 

patible with the experimental results of eq. (0). 

In parametrizing the matrix elements of the dimension-six current-current operators, 
we follow the analysis of ref. and we distinguish two cases, depending on whether or 
not the light quark q of the operator enters as a valence quark in the external hadronic 
state. Correspondingly, we introduce different 5-parameters for the valence and non- 
valence contributions. For the S-meson matrix elements, we write the matrix elements of 
the non-valence operators in the form 

{B,\0i\B,) _ flMs, 



1Mb 2 
while, in the case of the valence contributions (g = g'), we write 



br forg^g', (19) 



(B.mB,) flMs, {B,\Ol\B,) flMs, 



(20) 



The parameters 5^'^ in eq. (|20|) are defined as the 5^'^ of eq. (|T^) in the limit of degenerate 
quark masses (m^ = rriqi). In the VIA, Bi = B2 = 1 while the e parameters and all the 
5s vanish. Note that, in the SU{2) limit, the parameters B^2 express the matrix 

elements of the non-singlet operator (9^ — Of between external S-meson states. 

The reason to distinguish between valence and non valence contributions, is that only 



the former have been computed so far by using lattice QCD simulations 0, 10|. A non 



perturbative lattice calculation of the 6 parameters would be also possible, in principle. 
However, it requires to deal with the difficult problem of subtractions of power-divergences, 
which has prevented so far the calculation of the corresponding diagrams. 

To complete the definitions of the i?-parameters for the 5-mesons, we introduce a 
parameter for the matrix element of the penguin operator 

{B,\Op\Bl^j%Ms^^_ ^21) 



2M 



Ba 



We now define the 5-parameters for the A;, baryon. Up to l/rrib corrections, the matrix 
elements of the operators O2 and Ol, between external A;, states, can be related to the 
matrix elements of the operators Of and 

{Ab\Ol\A,) = -2 {AtmA,) , {A,\Ol\A,) = -2 (A,|0||A,) . (22) 

For the independent matrix elements, assuming SU{2) symmetry, we define 

2Ma, 2 V ^ ^ ^ 

2M.. = ^ + ^2 ) forg = t.,rf. 
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(OT.)^/|M,^,, for,^.,c, (23) 



O2 tor q 



2M, 



Ah 



(A,|0p|A,) flMs 
2Ma, 2 



P 



A 



In analogy with the i?-meson case, the parameters Li and L2 represent the valence contri- 
butions computed so far by current lattice calculations [|], |]. 

In terms of valence and non valence S-parameters, the quantities A^^g^ of eq. (p!^), 
which represent the spectator contributions to the lifetime ratios, are expressed in the 
form 



A^^ 

spec 


= 487r2 


PbMb 

^3^(3) 


spec 


= 487r2 


PbMb 


A^ 

spec 


= 487r2 


PbMb 

^3^(3) 



>k > 



k=l 

E [r 4 - + {cJ: + ct) (r ^ - 5^) + 

k=l 

{r 5r - St) + ct (r 5r - ^r)] +cp[rP^~P')] , (24) 

E [[4 + c.') ^fc^ - ct Bt + (c," + ct) (dt" - 5t) + 
fc=i 

(5^^^ - 5f ) + 4 {5t - 6^)] +cp{P^-P')}. 



where r denotes the ratio {f'^^^Ms^) / {f%MB) and, in order to simplify the notation, we 
have defined the vectors of parameters 

& = {BlBlelel}, 

£ = {Li,-Li/2,L2,-L2/2}, (25) 
5A^ = {5f'',-5f72, 52^-52^72}. 

An important consequence of eq. ( P^ ) is that, because of the SU{2) symmetry, the non- 
valence {5s) and penguin (Ps) contributions cancel out in the expressions of the lifetime 
ratio T{B~^)/T{Bd). Thus, the theoretical prediction of this ratio is at present the most 
accurate, since it depends only on the non-perturbative parameters actually computed by 
current lattice calculations. The prediction of the ratio r(A6)/r(Prf), instead, is affected 
by both the uncertainties on the values of the 6 and P parameters, and by the unknown 
expressions of the Wilson coefficients cj: and cp at the NLO. For the ratio t{Bs) /r^Bd) 
the same uncertainties exist, although their effect is expected to be smaller, since the 
contributions of non-valence and penguin operators cancel, in this case, in the limit of 
exact SU{?>) symmetry. 

In the numerical analysis of the ratios t{Bs) /r^Bd) and T{Ab) /r^Bd), we will neglect 
the non-valence and penguin contributions. The non-valence contributions vanish in the 
VIA, and present phenomenological estimates indicate that the corresponding matrix el- 
ements are suppressed, with respect to the valence contributions, by at least one order 
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of magnitude PB| , P7[ . On the other hand, the matrix elements of the penguin operators 
are not expected to be smaller than those of the valence operators. Since the coefficient 
function cp vanishes at the LO, this contribution is expected to have the size of a typical 
NLO corrections. Thus, from a theoretical point of view, a quantitative evaluation of the 
non-valence and penguin operator matrix elements would be of the greatest interest to 
improve the determination of the lifetime. 

By neglecting the non valence and penguin contributions, we obtain from eq. (0) the 
following numerical expressions for the spectator effects 

Afpt, = - 0.06(2) Bf - 0.010(3) + 0.7(2) ef - 0.18(5) , 

Ag%^ = - 0.010(2) Bf + 0.011(3) Bl - 0.16(4) + 0.18(5) e| 

+ 0.008(2) Bf - 0.008(2) B^ + 0.16(4) ef - 0.16(4) , (26) 

^^pec = - 0.08(2) Li + 0.33(8)^2 

+ 0.008(2) Bf - 0.008(2) B^ + 0.16(4) ef - 0.16(4) , 

These formulae, which are accurate at the NLO, represent the main result of this paper. 
The errors on the coefficients take into account both the residual NNLO dependence on 
the renormalization scale of the AB = 1 operators and the theoretical uncertainties on 
the input parameters. To estimate the former, the scale /ii has been varied in the interval 
between mi,/2 and 2mf,. For the charm and bottom quark masses, and the B meson decay 
constants we have used the central values and errors given in table ^ The strong coupling 
constant has been fixed at the value ^^(m^) = 0.118. The parameter c^^-* in eq. (|2^) is 
a function of the ratio ml /ml, and such a dependence has been consistently taken into 
account in the numerical analysis and in the estimates of the errors. For the range of 
masses given in table 0, c^^-* varies in the interval c'-^^ = 3.4 4.2 p2|-||24||. 

As discussed in the previous section, the HQE for the ratio t{B^) /T{Bd) can be also 
expressed in terms of operators defined in QCD. The corresponding coefficient functions 
can be evaluated by applying the matching defined in eq. ([T5|). In this way, we obtain the 
expression 

Afpt, = - 0.05(1) Bf - 0.007(2) B^ + 0.7(2) e'i - 0.15(4) 4 (27) 

where the B and e parameters are now defined in terms of matrix elements of QCD oper- 
ators. This expression is in agreement with the result obtained in ref. [^. 

The errors quoted on the coefficients in eq. (^) are strongly correlated, since they 
originate from the theoretical uncertainties on the same set of input parameters. For this 
reason, in order to evaluate the lifetime ratios, we have not used directly eq. (^61). Instead, 



we have performed a bayesian statistical analysis by implementing a short Monte Carlo 
calculation. The input parameters have been extracted with fiat distributions, assuming 
as central values and standard deviations the values given in table |^. The results for the 
5-parameters are based on the lattice determinations of refs. 0-[|l9-Q discussed in 



details in ref. P| , we have included in the errors an estimate of the uncertainties not taken 
''For recent estimates of these matrix elements based on QCD sum rules, see refs. p8|-[^. 



13 



Bf 


= 1.2 ± 0.2 


Bf-- 


= 1.0 ± 0.2 


Bl 


= 0.9 ± 0.1 




= 0.8 ± 0.1 


ei = 


= 0.04 ± 0.01 


= 


0.03 ± 0.01 


4 = 


= 0.04 ± 0.01 


^2 = 


0.03 ± 0.01 


Li - 


= -0.2 ±0.1 


L2 -- 


= 0.2 ±0.1 


rrib = 


4.8 ±0.1 GeV 


mb-rric = 


3.40 ± 0.06 GeV 


fB = 


200 ± 25 MeV 


fss / fs 


= 1.16 ±0.04 



Table 2: Central values and standard deviations of the input parameters used in the nu- 
merical analysis. The values of mi, and rric refer to the pole mass definitions of these 
quantities. 



into account in the original papers. The QCD results for the B meson i?-parameters of 
ref. |]10[ have been converted to HQET by using eq. (p!5|).f| The contributions of all the 6 



and P parameters have been neglected. 

In this way we obtain the NLO predictions for the lifetimes ratios which have been also 
quoted in the introduction 

^^^^} = 1.06 ±0.02, 4^ = 1-00 ±0.01, 4^ = 0.90 ± 0.05 . (28) 
r{Bd) T{Bd) T[Bd) 

The central values and errors correspond to the average and the standard deviation of the 
theoretical distributions. These distributions are shown in fig. ^, together with the experi- 
mental ones. We mention that uncertainties coming from the residual scale dependence in 
the results of eq. ( p8l) represent less than 20% of the quoted errors. 

In conclusion we find that, with the inclusion of the NLO corrections, the theoretical 
prediction for the ratio t{B^) /r^Bd) turns out to be in very good agreement with the ex- 
perimental measurement, given in eq. (0). For the ratios t{Bs) /T{Bd) and T{Ah)/T{Bd) the 
agreement is also very satisfactory, and the difference between theoretical and experimen- 
tal determinations is at the la level. We have pointed out, however, that the theoretical 
predictions are less accurate in these cases, since a reliable estimate of the contribution of 
the non- valence and penguin operators cannot be performed yet. We also found that the 
NLO charm quark mass corrections computed in this paper are rather large for some of 
the Wilson coefficients. Nevertheless, the total effect of these contributions on the lifetime 
ratios is numerically small. 



^With respect to ref. we use for the i?- meson _B-parameters the results updated in ref. ||lC 
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Figure 4: Theoretical (histogram) vs experimental (solid line) distributions of lifetime ra- 
tios. The theoretical predictions are shown at the LO (left) and NLO (right). 
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Appendix 

In this appendix we collect the analytical expressions of the Wilson coefficients, both at 



the LO and NLO. The LO coefficients have been computed in refs. and are reported 

here for completeness. 

We distinguish the leading and next-to-leading contributions in the coefficients F^^j by 
writing the expansion 

= + ^Bl^ ' (29) 

(g = u,d',s'). Since, by definition, the coefficients Al ^j and Bl ^j are symmetric in the 
indices i and j, we will only present results for i < j. 



The LO coefficients Al read 



^1.11 - 3 

^2,11 — ' 



A- = 2 (1 - z)^ 



4" — n 



_ -i^-zf (2 + z) 



Ad' _ ~j y- I ~, 

^1,11 — 2 ' ^1,12 

= (1 + 2^) , <12 



Ad' _ n Ad' 
^3,11 ~ ' ^3,12 



^4,11 — ! ^ 



4,12 



AU 

^1,12 


^{l-zf, 


All _ 
^1,22 — 


a-z] 
3 


7 




AU 

^2,12 


= 0, 


All _ 
^2,22 ^ 


0, 






AU 

^3,12 


= 0, 


AU _ 
^3,22 — 


2(l-zf , 




AU 

^4,12 


= 0, 


AU _ 
^4,22 — 


0. 






_ -(1 


-zf (2 + z) 


Ad' 
^1,22 


_ -(1 


-zf (2 


+ z) 




6 




18 




_ (1- 


zf {l + 2z) 


Ad' 
^2,22 


(1- 


zf (1 + 


2z) 




3 




9 




= 0, 




Ad' 
^3,22 


_ -(1 


-zf (2 


+ z) 






3 




= 0, 




Ad' 
^4,22 


2(1 


-zf (1 


-f 2z) 






3 





(30) 



(31) 



16 



^1,11 
Ai 



-VI -4z {1-z) 



^2,11 



^3,11 



VI -4z (l + 2z) 

0, 

0, 



1-2,12 



^1,12 

A1 



-VI -4z 
3 

Vl-4z (l + 2z) 
3 

0, 
0, 



^1,22 
^2,22 
A?' 



-VI -4z (1-z) 
9 

Vl-4z (l + 2z) 
9 



-2Vl-4z (1 - z) 



2 Vl - 4z (1 + 2z) 



where z = m^/m^. 



(32) 



The NLO results for the coefficients - have been obtained in the NDR-MS scheme 



of ref. [0] for the AS = 1 operators and the NDR-MS scheme of ref. for the Ai? = 0, 
HQET operators. We find 
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-8 (1-z) (164 + 127r2 - 109z + 5z2) 16 (1 - z)^ (6 logxi + (1 + z) log(l - z)) 



81 9 

16z (l2 + z-3z2) logz 16 (l-3z) (1-z) (log(l - z) log z + 2 Li2(z)) 



27 



9 



B. 



241 



-^^3,11 



16 (1 - z) (4- llz+ 10z2) 16z2logz 

81 ^ 27 ' 

((1 - z) (751 - 127r2 (7- 9z) - 1199z- 26^2)) 

27 

,2 4 (26 - z) (1 - zf log(l - z) 
8 1oga;i(l-z) 

2z (192- 119z + 6z2) logz 8 (4-3z) (1 -z) (log(l - z) log z + 2 Li2(z)) 

9 3 
-8 (1-z) (ll-10z-13z2) 32z2logz 



27 



(33) 



-8 (1-z) (l3 + 27r2-2z + z2) 4(l-z)^ (6 logxo + 4 log(l - z)) 

9 3 

16 (3-z) zlogz 16 (1-z) z (log(l-z) logz + 2Li2(z)) 



Si 



-°3,12 



^4,12 



32(1 - zy 
9 ' 

-2 (1 - z) (125 + 6 TT^ - 103 z + 2 z^) 

9 ^ 

x2 . . ^^ 4z (3 + 4z-3z2) log; 
(1 - zf (6 loga;o - 24 logxi - 4z log(l - z)) ^ '—^ 

4 (1-2 z) (1-z) (log(l-z) logz + 2Li2(z)) , 
-4 (1 -z) (2-z-4z2) 4z2logz 



(34) 



~8 (1 - z) (164 + 12 - 109 z + 5 ; 



81 
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16 (1 - zf (6 log xi + (1 + z) log(l -z)) 16 z (12 + z - 3 log 2 

9 27 
16 (1-3^) {1-z) (log(l--g) log-2 + 2Li2(-z)) 
9 

16 (1-z) (4-11 z+ 10^2) 16^2 logz 

81 27 ' 

(1-^;) (751 -875z-26 02_i2 7r2 (7 + 9^)) 
27 

4(1- z)^ (6 logxi - (17- z) log(l - z)) _2 (28- 3 z) (3-2z) z logz 

3 9 
8 (1-z) (5 + 3 z) (log(l-z) logz + 2Li2(z)) 
3 

-8 (1-z) (ll-10z-13z2) 32z2logz 
27 9 ' 



2 (1-z) (2-z) (5 + 3z) 4(l-z)^ (4 + 5z) log(l - z) 

3 ^ 3 ^ 

4z (5 + 4z - 5z2) logz 



+ 

6 

2 (1 - zf (2 + z) (6 loga;o - 4 (log(l - z) logz + 2Li2(z))) 

3 

-4 (1-z) (5-7z + 6z2) 8(l-z)^ (2 + 10z-3z2) log(l - z) _^ 



3 

8z (2 - 17z + 16z2 - 3z3) logz 



4 (1-z)' (l + 2z) (6 loga:o-4 (log(l-z) log z + 2 Li2(z))) 

3 

-((1-z) (205 + 7Z- llOz^)) 



18 

3(1 - zf (2 + z) (2 logxo - 2 log(l -z)) _ z^ (6 + Hz) logz 

2 3 ' 

(1-z) (95 + 104^-211z2) 
9 

(1 - zf (1 + 2 z) (6 log.o - 6 log(l - z)) - 4(12-ll ^)^Mogz 



!7r2 (1 - z) z (2 + z) ^ 2 (1-z) (42 - 15 z - 25 z^) ^ 



27 

2 (1 - zf (2 + z) (2 logXQ + 4 logxi) 
3 

4(l-z)' (2 + z + 3z2) log(l-z) 4z (l + 6z-6z2) logz 

9z ^ 9 

8 (1 - z) (2 - z) (2 + z) (log(l - z) logz + 2Li2(z)) 

9 

-167r2 (1-z) z (l + 2z) 4 (1-z) (21 + 27z - 38 z^) 



27 9 
4 (1-z)' (l + 2z) (2 1oga;o+4 1oga;i) _ 
3 
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8{l-zf {l + 2z + 6z^ -3z^) log(l - z) 
9^ 

8z (4 - 24 2; + 18 0^ - 3 2;^) logz 

9 ^ 
16 (1 - z) (2 - 2) (1 + 2 z) (log(l -z) logz + 2 U2{z)) 

9 

2 7r^ (1-2) ^ (2 + ^) (1-2) (83-1512-88^^) _ 

9 54 
(1 - zf (2 + z) (2 logXQ - 4 logxi) 
2 

(1 - zf (1 + 2) (2 + z) log(l - z) _ 2 2 (6 + 72^) log2 

32 9 
2 (1 - 2) (2 + 2) (log(l - 2) log 2 + 2 Li2(2)) 

3 

-47r2 (1-2) 2 (I + 22) (1-2) (49 + 202 2- 185 2^) 



+ 



9 27 

(1-2)^ (1 + 2 2) (2 loga;o-4 1oga;i)- 

2(1- zf (1 + 2) (1 + 2 2) log(l - 2) 22 (6-452 + 28 2^) log2 

3~z 9 
4 (1-z) (I + 22) (log(l-2) l0g2 + 2Li2(2)) 



167r' jl-z) z (2 + 2) 2 (1-z) (461 -286 2-313 2^) 

81 243 ^ 

16(1-2)^ (2 + 2) logai 16(1-2)^ (1 + 2 + 2^) log(l-2) 
9 272 

4 z (9 - 18 2 + 32 2^) log 2 
81 

8 (1-2) (3-2) (2 + 2) (l0g(l-2) l0g2 + 2Li2(2)) 

27 

-32 TT^ jl-z) z (1 + 22) 4 (1-2) (238 + 445 2-5272^) 

81 243 
32(1-2)^ (1 + 22) loga;i _ 
9 

8(1-2)^ (2 + 52 + 82^-323) log(l-2) 

27^ 

82 (18- 1172 + 82 22-923) log2 



+ 



+ 



81 

16 (1-2) (3-2) (1 + 2 2) (l0g(l-2) l0g2 + 2Li2(2)) 

27 

27r2 (1-2) (2 + 2) (9 + 72) (1-2) (937 -449 2- 746 2^) 
27 81 

4(1-:)' (2 + c)logxi 4(l-r)- (1-17--8--) los(l - 2) 

3 92 
2z (72 - 92- 20 2^) log2 
27 

4 (1 - 2) (2 + 2) (3 + 5 2) (log(l - 2) log 2 + 2 Li2(2)) 
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-47r^ {1-z) {l + 2z) {9 + 7z) {1 - z) (715 + 1003 z - 2804 z^) 

27 81 
8(1-^)' (1 + 2^) loga;i 
3 

2(l-z)^ (2-31^-64^2-323) log(l-z) 
9^ 

2z (36-288^ + 62 02^9^3) jQg^, 

27 ^ 
8 (1 - z) (1 + 2 z) (3 + 5 z) (log(l - z) log z + 2 Li2(z)) 



4V1-40 jl-z) (5 + 6-g) 8 (6- 13 2:-202^g_g3^ igg^ 

4V1 -4^ (1 - z) (6 logxo + 8 log(l - 4z) - 12 logz) 

3 

^ (1 - 2 2;) (1 - 2;) (3 logV + 2 log cr log(l -4^) - 3 log cr log0 + 4Li2(CT)+ 

4V1-40 (1-4-2) (5 + 6-g) _ 16 {3-2z-7z'^ + 12z^) logo- ^ 
16(1-2 .) (1 + 2 .) log^ . + ''^'-'^^ (l + 2!)logolog(l-4.) _ 



4V1-40 (1 + 2^) (6 1oga;o + 8 1og(l-40)-12 log-g) _ 

16 (1-2Z) (1 + 2.) logoiL + ^^iii^M(i±MiM^ + 

o 

32 (1-22) (1 + 2 0) Li2(o2) 
3 ' 

- (VI -42 (205 + 14 2 + 24 22)) 2 (9 - 272 - 622 + 42^) logo- 



18 3 
3\/l - 4z (1 - z) (2 logxo - 2 log(l - 4z) + 2 logz) , 
Vl-42 (95 + 208 2 + 48 22) 2 (9-622 + 162^) logo 

3\/l -42 (1 + 22) (2 1oga;o-2 1og(l-4 2) + 2 1og2) , 



\/l-42 (26- 172 - 622) (1 -602+14622 - 362^) logo- _^ 



+ 



3 9z 
4 VI -42 (1 - 2) (6 logxo + 12 logxi +8 log(l - 42)) 

9 

Vl-42 (1 - 582 + 72 22) log2 
92 

y (1 - 22) (1 - 2) (3 log^o- + 2 logo- log(l -42) - 3 logo- log2 + 4Li2(o)+ 
2Li2(a2)) , 

-4 Vl-42 (7-32) (1 + 22) _^ 4 (1- 152-422 +4823- 362^^) logo- 



3 92 
4V1-42 (1 + 22) (6 loga;o + 12 loga;i + 8 log(l - 4 2)) 

9 

20 



4 VI -42 (1 - 13 2- 36^2) log^ _^ 



9^ 

y (1 - 2 (1 + 2 2) (3 logger + 2 logo- log(l - 4z) - 3 logcr log^ + 4Li2((T)- 

2Li2((72)) , 

y/l-Az (ll2-5232 + 6z2) (3-108^ + 342^2^4 04) log cr 
108 36z 



Vl - 4z (1 - z) (2 log .To - 4 logxi - 2 log(l - 4^)) + 
Vl-4z (1-340 + 48 2^) logz 

120 ' 
- (VI -40 (49 + 2150 + 6 0^)) (3- 270-3602 + 7203 + 40^) loga 



27 9 



+ 



Vl - 40 (1 + 2z) (2 logxo - 4 loga;i - 2 log(l -4 0)) 
Vl-40 (1 - 70-2402) logz 



30 



2 VI -40 (407 - 4910 - 78 02) 2 (3 - 144 + 390 02 - 52 0^) logcr 

243 81^ ^ 

8V1 -40 (1-0) (12 1oga;i+7 1og(l-4 0)) 2 VI -40 (1 - 460 + 6002) log; 



+ 



27 27^ 
^ (1-20) (1-0) (3 1ogV + 2 1og(7log(l-40)-3 1ogCTlog0 + 4Li2((T)+ 
2Li2(a2)) , 

-8 Vl-40 (119 + 256 - 78 02) 8 (3 - 360 - 2402 + 1080^-5204) loga 

243 ^ 81^ 

8 Vl-40 (1 + 20) (12 1oga;i + 7 1og(l-4 0)) _ 8V1-40 (l- 100 - 30 02) log0 
27 270 

^ (1-20) (1 + 20) (3 log2a + 2 logtT log(l-4 0)-3 logtT log0 + 4Li2(a)+ 
2Li2(a2)) , 

47r2(l-0) Vl-40 (1129 - 2143 + 186 z^) 

3 81 
2 (21 + 780- 267 02 + 90 0^ + 62 04) log cr 
27^ 

8 (2-50) Vl-4^;loga;i 112 VI -40 (1 - z) log(l-40) _ 

9 9 
4 (1-0) (7 + 40) log cr log(l -40) _ 2 Vl-40 (7 + 400- 7102) log0 

9 90 
4 (1-0) (5 + 20) (log2 cr- logo- logz) ^ 
3 

16 (1 - 2 0) (1 - 0) Li2(o-) _ 16 (1 - 0) (4 + 0) 112(0-2) 

9 9 ' 

-4 7r2 (1 + 2 0) Vl-40 (691 + 1922 + 1608 02) 

3 ^ 81 

2 (3 + 159 + 138 02 - 774 + 536 0^) log a 



27 

8V1 -40 (1 + 20) (2 loga;i-14 log(l-40)) 



+ 

+ 
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4 (1 + 2z) (7 + 4z) logcr log(l -4z) 2^/1 -4z (l + 55 z + 154 z^) logz 

9 ^ 9~z ^ 

4 (l + 2z) (5 + 2z) (log2 cr - log cr logz) 

(41) 



3 

16 (1 - 2 z) (1 + 2 z) Li2(cr) 16 (4 + z) (1 + 2 z) Li2(cr2) 



9 9' 
where a is the ratio 

g^^-^, (42) 

and we have defined = fio/^b and xi = fii/nib. 

Finally we present the results for the coefficients Pk^ij of the penguin and chromomag- 
netic operators defined in eq. ([l^) . The coefficients Pk,2s have been computed by using the 
convention in which the chromomagnetic coefficient CsG has a positive sign. We obtain 
the expressions: 



\/r^4i(i-z) „ -yr^4i(i-z) 



Pl.l3 = -VT^4i (1 - z) , Pi,23 - ^ , Pl 



Vl-4z (l + 2z) Vl-42 (l + 2z) 



P2.13 = \/r^^(l + 2z) , P2 23 ^ , P2,14= ' ^ (43) 

^^3,13 = , P3.23 = , P3,14 = , 

^4,13 = , P4,23 - , P4,14 = . 



_ -^/^^4I (1-z) 

-fl,24 — 7: 



_ VT^4l (l + 2z) 

-f2,24 — T 



_ -2VT^4l (1-z) 

-f3,24 — 7, 



2V1 -4z (l + 2z) 



A, 15 


= -3zVl-4z, 


Pl.25 


= -z Vl - 4z 


^2,15 


= 0, 


-^2,25 


= 0, 


-P3,15 


= 0, 


-^3,25 


= 0, 


^4,15 


= 0, 


^4,25 


= 0. 



(44) 



Pl,16 = -Z\/r^4l, Pl,26 = ^^^L^ , Pl^28 = 0, 

^2.16 = , P2.26 = , P2,28 — , 

P3a= = 0, P3,» = -2 = ^, _-2^/mI(l + 2,) (45, 



3,28 



2V1 -4z (l + 2z) 

^^4,16 - , P4,26 = , P4,28 = — ' 



Note that, in the limit of vanishing charm quark mass {z = 0), the contribution of the 
penguin operators and Qq vanish for chirality. 
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